ASYMPTOTIC BEHAVIOR FOR SYSTEMS OF NONLINEAR 
WAVE EQUATIONS WITH MULTIPLE PROPAGATION 
SPEEDS IN THREE SPACE DIMENSIONS 



SOICHIRO KATAYAMA 

Abstract. We consider the Cauchy problem for systems of nonlinear wave 
equations with multiple propagation speeds in three space dimensions. Un- 
der the null condition for such systems, the global existence of small am- 
plitude solutions is known. In this paper, we will show that the global 
solution is asymptotically free in the energy sense, by obtaining the asymp- 
totic pointwise behavior of the derivatives of the solution. Nonetheless we 
can also show that the pointwise behavior of the solution itself may be quite 
different from that of the free solution. In connection with the above results, 
a theorem is also developed to characterize asymptotically free solutions for 
wave equations in arbitrary space dimensions. 



I. Introduction 

We consider the Cauchy problem for a system of nonlinear wave equations 
of the following type with small initial data: 

D Cj Uj(t,x) = F j (du(t,x),d 2 u(t,x)), (t,x) e (0, oo) x IR 3 , (1.1) 

Uj(0, x) = efj(x), (d t Uj)(0,x) = egj(x), 16K 3 (1.2) 

for j = 1, . . . , N, where D c = d 2 — c 2 A x = d 2 — c 2 J2k=i °l k f° r c > 0, and 
the propagation speeds ci, . . . , cjy are positive constants, while du and d 2 u 
denote the first and second derivatives of u = (w/)i<Kjv, respectively. More 
specifically, we write 

du = (d a Ui)i<[<N,0<a<3, 9 2 U = (<9 a <9feWi)i<K7v,0<a,b<3 

with the notation 

d d 

d := d t = and d k := d Xk = —— for k = 1, 2, 3. 
at ox k 

We assume that / = (fj)i<j<N,9 = (9j)i<j<N e C °°(M 3 ; R N ). e in (OJ) is a 
small and positive parameter. Fj can include nonlinear terms of higher order, 
but for simplicity we assume that Fj can be written as 

N / 3 3 \ 

F > = E E Pfi (daU k )(d b d b , Ul ) + ^ 4l(daU k )(d bUl ) (1.3) 
k,l=l \o,6,6'=0 a,b=0 / 
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for 1 < j < N with appropriate real constants and qfy, where p°j$ = 
for 1 < j,k,l < N and < a < 3. To assure the hyperbolicity, we assume the 
symmetry condition 

P%=pt=pfi (1-4) 
for 1 < j, k, I < N, and < a, b, V < 3. 

For general quadratic nonlinearity, the classical solution to f ll.ip -f fL2|) may 
blow up in finite time for some (/, g) no matter how small e is. Klainerman 
[T6] introduced a sufficient condition for small data global existence for the 
single speed case where C\ — • ■ ■ — c^(= 1) (see also Christodoulou [2]). This 
sufficient condition, called the null condition, was extended to the multiple 
speed case: To simplify the description, we assume that the speeds are distinct, 
and that 

< ci < c 2 < • • • < Cat. (1.5) 

Let the constants and gjjy be from (11. 3p . We say that the null condition 
(associated with the speeds ci, . . . , cn) is satisfied if we have 

3 3 

pf-X a X b X bl = '/;.!'rV,.V,, = 0, X G A/}, 1 < j < N, (1.6) 

a,b,b'=0 a,b=0 

where Mj := {X = {X a ) < a <3 G M 4 ; X 2 - Cj 2 ^Li X k = °>- Sma11 data S lobal 
existence under the null condition for the multiple speed case was obtained 
by Yokoyama [29] (see also Sideris-Tu [27], Sogge [28], and Hidano [1]; see 
Kubota-Yokoyama [TS], the author [7J, and Metcalfe-Nakamura-Sogge [22] for 
the case where nonlinearity of higher order depends not only on (du, d 2 u), but 
also on u). We introduce the null forms 

3 

QoM;c) =(dtV)(d t ip) -c 2 Y,(dk^(d^), (1.7) 

fc=i 

Q ab ( V ^)=(d aV )(d b ij)-(d bV )(d a ^), 0<a,6<3. (1.8) 

Then it is shown in [29] that the quadratic nonlinearity satisfying the null 
condition can be written as 

Fj(du, d 2 u) = Nj(du, d 2 u) + R)(du, d 2 u) + Rf(du, d 2 u), (1.9) 

where 

3 

N i= Yl ^Q {u h d a uf lCj )+ Bf b 'Q ab (u 3 ,d b * Uj ), (1.10) 

0<|o|<l a,b,b'=0 

R )= E ( E Pfki (d a u k )(d b d b/Ul ) + ]T qf kl (d a u k )(d bUl ) ] , (1.11) 
R f= E ( E pS'(W(^ Wfe )+^g^ fc (^)(^)] (1.12) 
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with some constants A? and Bf bb ' (p"|f and q°- h kl are from (11. 3 ft ). Here d = 
(d a )o< a <3, and a = (a a )o<a<3 is a multi-index. We refer to terms involved in 
R 1 = (R))i<j<N and R 11 = (Rj)i<j<N as the non-resonant terms of type / and 
type II, respectively. 

Now we turn our attention to the asymptotic behavior of the global solutions. 
Let c > 0. It is known that if G G ^ 1 ((0, oo); L 2 (M 3 )), then the solution v to 
D c v(t,x) = G(t,x) is asymptotically free in the energy sense, that is to say, 
there exists a solution v + to the free wave equation U c v + = such that 

lim \\(v-v + ){t)\\ E ,c = 0, 

t— >oo 

where the energy norm ||<£>(i)||E,c is defined by 

Mt)\\ 2 E,c = l^ 3 t^\dMt,x)\ 2 + \VMt,x)\ 2 ^j dx 

with V x = (di,d2,d 3 ). For the single speed case, investigating the proof in 
[TB], we have 

Fj(du, d 2 u) G L 1 ((0, oo); L 2 (R 3 )) 

under the null condition because of the extra decay for the null forms, and as 
an immediate consequence we see that the global solution u is asymptotically 
free in the energy sense. As for the multiple speed case, by the estimates 
obtained in [18] (cf. (15. 6p and (15.7ft below), it is easy to see that 

N j (du,d 2 u) + R)(du,d 2 u) G L 1 ((0, oo); L 2 (IR 3 )) (1.13) 

for the global solution u to (ll.ip -( !TT2|) . because we can expect some gain in the 
decay rate for the null forms and the non-resonant terms of type I, compared 
to general quadratic nonlinearity (see (I5.13P and (I5.14p below). Therefore, if 
R n = in (II. 9p . then (j 1.13p implies that the global solution u is asymptotically 
free in the energy sense; namely for j = I, . . . , N, there exists a solution 
of the free wave equation O c .u^ = such that lim^oo \\(uj — u^)(t)\\ E)C:j = 0. 
By contrast, there is no explicit gain in the decay rate for the non- resonant 
terms of type II unless they can be written in terms of the null forms, and we 
cannot expect R^(du, d 2 u) G L 1 ((0, oo); L 2 (M 3 )) in general (see (15.15ft below), 
although its influence is weak enough for the solution to exist globally. Hence 
it is not clear whether the global solution u is asymptotically free or not when 
R 11 is not written in terms of the null forms. Our aim in this paper is to 
determine the asymptotic behavior in the presence of the non-resonant terms 
of type II by modifying the method developed in [10J. 



2. Main Results 

2.1. Asymptotically free functions in the energy sense. To begin with, 
we will characterize the asymptotically free functions in the energy sense. We 
only need the three space dimensional result for an application in this paper, 
but we consider the general space dimensional case here for future applications. 
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Let n be a positive integer. Let i7 1 (R n ) be the completion of C^R") with 
respect to the norm H^lljjirRn) = II ^x<p\\L 2 (R n )- We put 

X n := ^([O.oo);^ 1 ^™)) fl C 1 ([0, oo); L 2 (R™)). (2.1) 

We say that a function v = v(t, x) G X n is asymptotically free in the energy 
sense associated with the speed c, if there is (vq,Vi) G iy 1 (R n ) x L 2 (R") such 
that 

lim \\v(t,-)-v + (t,-)\\ E , c = 0, 

t—KX) 

where v + G X n is a unique solution to 

n c v + (t,x) = 0, (i, x) G (0, oo) x R n 

with initial data (t> + (0), d t v + (0)) = (vq, vf), and the energy norm (associated 
with the speed c) is given by 

h(t, OIII,c : = ^ £ *)l 2 + |v,v(t, *)l 2 ) cfa. 

Here D c = <9 2 — c 2 Ylk=i ®1 with the notation d k = d/dx k for 1 < k < n. Note 
that we do not suppose that v is a solution to some wave equation here. 
For c > and x G R n \ {0}, we define 

U c (x) = iyJc,a( X )) < a < n : = (~C, M^)- (2.2) 

Theorem 2.1. Lei n > 2 and c > 0. ^4 function v G X n zs asymptotically 
free in the energy sense associated with the speed c, if and only if there is a 
function V = V(cr,u) G L 2 (R x S"- 1 ) such that 

&IN*.o -^(0^,011^ = 0, 

where d = (d , d\, . . . , d n ), and 

V\t, x) := \x\- {n - 1)l2 V{\x\ - ct, Ixl^x), (t, x) G (0, oo) x (R n \ {0}). 

This theorem for n = 3 was implicitly proved in (TO]. The author believes 
that this result for general space dimensions is a new observation, but it possi- 
bly has already appeared in some literature. We will give the proof in Section[3j 
See [H] for an application of Theorem 12. II to semilinear wave equations in two 
space dimensions. Another application in three space dimensions can be found 
in [13]. 

2.2. Asymptotic behavior of the global solutions. Next we examine the 
asymptotic behavior of the global solutions to fll.ip - fll.2l) . For h G C°°(R 3 ), 
we define its Radon transform lZ[h] by 

n[h}(a,u)= [ h(y)dS(y), 

where dS(y) denotes the surface element on the plane {y ■ u = a}. Now, 
restricting our attention to the three space dimensional case, we introduce 
the Friedlander radiation field (its definition in general space dimensions will 
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be given in (ET25]) below): For ((p,ip) G C£°(R 3 ) x C£°(R 3 ), we define the 
Friedlander radiation field 

TbfoVK*,") = ^{n[iP}(a,u) - (d a K&])(a,u)) (2.3) 

for (o~,ou) G R x § 2 . For z G R d with a positive integer <i, we put (z) := 
y/T+W- 

Theorem 2.2. Suppose that (jl.3p . (11.4)1 . ( II .5|) . and iae nit// condition ( II. 6p 
are fulfilled. Let d = (d , di, d 2 , d 3 ), and cu c (x) be given by ( 12.2)1 . 

(1) We sraali o" > 0. Tnen /or any /, a G C^°(R 3 ; R^) and sufficiently small 
e > 0, there is a function P = (Pj)i<j<N of (a,u) G R x § 2 such that 

\x\duj(t,x) =eco Cj (x)Pj(\x\ — Cjt, |x| _1 x) 

+ 0(s(t+\x\)" 1+s ( Cj t-\x\)- 5 ) (2.4) 

for (t,x) G [0, oo) x (R 3 \ {0}) and 1 < j < N, where u = (uj)i<j<N is the 
global solution to the Cauchy problem ( 11. ip - 111. 2p . Moreover we have 

Pj{a,u) = d a T [f j ,cJ 1 g j ){a,u;)+O(e{a)- 1 ), (a, to) G R x § 2 . (2.5) 

(2) We further assume that R 11 = {R l ^)i<j<N has the null structure, that is to 
say 

R f= J2 f A« k Q (u k ,d a u k ;c k )+ Bf'Q ah {u k ,d v u k )\ 

{k;k^j} \0<|a|<l a,b,b'=0 J 

for 1 < j < N with some constants A" k and B a ^' . Then for any f,g£ 
C£°(M. 3 ,M. N ) and sufficiently small e > 0, there is a function U = (Uj)\<j<N 
of (a, oo) G R x § 2 such that 

\x\uj(t,x) =eUj(\x\ - Cjt, \x\~ 1 x) + 0(e (t + log(2 + t)), (2.6) 
\x\duj(t,x) —£UJ c ,(x)(d a Uj)(\x\ — Cjt, \x\~ x) 

+ 0(e(t+\x\y 1 (c j t-\x\)- 1 ) (2.7) 

for (t,x) G [0, oo) x (R 3 \ {0}) and 1 < j < N, where u = (uj)\<j<N is the 
global solution to the Cauchy problem (jl.ip - (jl.2p . Moreover we have 

d™Uj(a,u) = d™Mf3,c- l g] ]{a,uj) + 0{E{o)- l - m ), (a,u) G R x S 2 (2.8) 

for m = 0, 1 and 1 < j < N. 

By ( 12 .4p we see that the asymptotic pointwise behavior of duj is similar to 
that of the derivatives of the free solution (see Lemma [3.31 below). Combining 
(12 .4p and ( 12 .5p with Theorem 12.1) we see that the solution u is asymptotically 
free in the energy sense. 

Corollary 2.3. Suppose that (II. 3p . ( II. 4p . (I1.5p . and the null condition ( II. 6p 
are fulfilled. Then, for any f,g G C^°(R 3 ;R iV ) and sufficiently small e > 0, 
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there exist f + = (f]~)i<j< N e H^M. 3 ;^) and g + = (gf)i<j<N e L 2 (R 3 ; R N ) 
such that 

lim\\( Uj -u+)(t)\\ EtCj = 0, l<j<N, 

where u = {uj)\<j<N is the global solution to the Cauchy problem (11.11) - (11.21) . 
and ?x+ is the solution to CLit^ = with initial data {u\,d t u\) = (/ + ,gt) at 



t = 0. 

Theorem 12.21 and Corollary 12.31 will be proved in Section [SJ after giving some 
preliminaries in Section HI 

When R 11 has the null structure, we see from (I2.6P that the solution Uj itself 
also behaves similarly to the free solution. The next result shows that the lack 
of the estimate corresponding to (I2.6P in the presence of R without the null 
structure is inevitable. 

Theorem 2.4. Let < c\ < C2, and let A±,A 2 be real constants. We consider 
the Cauchy problem for 



D Cl ui = Ai(<9 t u 2 ) 2 , 
n c2 u 2 = A 2 (d t u 1 ) 2 



in (0, oo) x R 3 . (2.9) 



If At ^ 0, then there exist (f,g) G C^°(IR 3 ; R 2 ) x C^^jM 2 ), M > 0, T > 2, 
and C > such that 

C"V(l + £log(2 + t)) < \r Ul (t,x)\ < Ce(l + Elog(2 + t)) (2.10) 

for any (t,x) satisfying T < c x t < r(= \x\) < c\t + M, where u = (mi,m 2 ) is 
the global solution to (12.91) with initial data u = ef and d t u = eg at t = for 
sufficiently small e(> 0). 

This theorem will be proved in Section [6j From this result, we see that 
there is some loss in the decay rate of u\\ ( I2.10p says that Ui(t, x) decays like 
(1 +t) _1 log (2 + t) along the ray r = c\t + a for < o < M, while (12. 6p implies 
that Uj decays like (1 + t)" 1 and has the same decay property as the free 
solution along the ray r = Cjt + a if R 11 has the null structure. Accordingly we 
find that estimates like ( 12 .6p cannot hold in the presence of the non-resonant 
terms of type II without the null structure. 

To sum up the results, an interesting character of the non-resonant terms 
of type II is revealed: Their effect is weak enough for the solution to exist 
globally (Theorem 1.1 in [29]) and to be asymptotically free in the energy 
sense ( Corollary 12. 3p ; however it is strong enough to affect the decay rate of the 
solution u itself (Theorem 12 .4p . though that of du is not affected ( Theorem 12 .2 j) . 

Throughout this paper, various positive constants are denoted by the same 
letter C. Thus the actual value of C may change line by line. 

3. ASYMPTOTICS FOR HOMOGENEOUS WAVE EQUATIONS 

Our aim in this section is to prove Theorem 12.11 Most of the necessary 
materials for this purpose are rather standard, but we give the details and 
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proofs to make this section self-contained. In what follows, we use the formal 
expression of writing distributions as if they are functions. 

3.1. Solutions to the free wave equation. We put H (R n ) = if^R") x 
L 2 (R n ), and we define 



IK^VOIIjW") := 




iJ (IR n ) can also be understood as the completion of C^°(lR n ) x C~(lR n ) with 
respect to the norm 1 1 * 1 1 n$n ^ . Let c > 0, and consider the Cauchy problem 
for the free wave equation 

D c w{t,x) = 0, (t, x) G (0, oo) x IR n , (3.2) 

(w(0,x), (d t w)(0,x)) = (w {x),wi{x)), xeR n . (3.3) 

Let X n be defined by (EH). It is known that for (w , wx) G H Q {R n ), fl3T2|) -f l3~3|) 
admits a unique solution w G X n , and we have the conservation of the energy 

IHV)IU,c= \\w(0, 0IU,c= ll(wo,c _1 Wi)||i/ (Rn). 
For Re a > — 1, let \+ be defined by 

( s a 

x a + (s)-.= {WTTy s>0 ' 

[0, s < 0, 

where T denotes the Gamma function. Then we have x+( s ) = f° r 
Re a > —1. We can extend the definition of the distribution x+ to all a G C 
so that we have x+( s ) — (x+ +1 )'( s ) f° r an y a £ C. Note that x+ can have its 
singularity only at s = 0. Especially we have 

X ;*(a)=«5 ( *- 1) (a), keN, 

where 5 is the Dirac function and 6^ denotes its j-th derivative. For a positive 
integer m, we define 

*»<*.*) = ^r,xt m " 2 (t 2 - N 2 )- 

Then the solution w to (|53jl - (|33)» with (tu ,Wi) e (C£°(M n )) 2 can be written 

as 

w(t, x) = c^d^Enict, •) * w ) (x) + c" 1 (E n (ct, ■) * (3.4) 

where the convolution * is taken with respect to x- variable (see Hormander [51 
Section 6.2] for instance). 

From now on, we suppose that (w ,wi) G (L7 °°(M n )) , and that 
Wi(x) = for \x\ > M with some positive constant M. Since supp E n (t, ■) C 
{x; \x\ < £}, it follows from (13 .4p that 

w(t,x) = 0, \x\ >ct + M. (3.5) 

If n(> 3) is odd, then supp E n (t, ■) = {x; \x\ = t} and we also get 

w(t,x) = 0, \x\<ct-M, (3.6) 
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which is called the (strong) Huygens principle. If n is even, (I3.6P is not valid 
in general, but we have a faster decay away from the light cone ct = \x\: Let 
t/2 > \x\ and t > AM, say. Then t 2 - \x - y\ 2 > 7£ 2 /16 > C (t + \x\} 2 > 
for \y\ < M. Hence, observing that x^~~ n ^ 2 {s) = A n s^ l ^ n ^ 2 for s > with an 
appropriate constant A n , we get 

\d« x E n {t,x-y)\<C a (t+\x\y lal+{1 ~ n \ t/2 > max{|x|, 2M}, \y\ < M 

with a positive constant C a . Therefore (13 .4p leads to 

\d a w(t,x)\ < C a {t+ |x|> HQ|+(1 ~ n) , ct/2>\x\, (3.7) 

because (13.71) for 2M > ct/2 > \x\ is easily shown. 

Following the arguments in Hormander [6], Section 6.2], we will obtain a 
useful expression of E m (t) * tp for ip G Cg°(R n ) with tp(x) = for |x| > M. 
Note that we have 

(E m {t) * ip) (x) = 0, \x\>t + M, (3.8) 

and 

(E m {t) * ip) (x) =0, \x\<t-M when m(> 3) is odd (3.9) 
as in (I3.5P and (13.61) . Let x = ruj with r = \x\ and u G S n_1 . We assume 

t 

2M < - < r < t + M. 
~ 2 ~ ~ 

We put a = r — t. Then we get — r < a < M. Since we have 

t 2 — \x — y\ 2 = 2r(u ■ y — a) + o 2 — |y| 2 , 
the homogeneity of x+~ m ^ 2 implies that 

(2nr)^/ 2 (E m (t) * ip)(x) =\ £ xt^ ■ y - a + <p(v)<fo 

where (?[<£>] (s, w, z) is given by 

Q[<p](,s,u,z) = J 5 (^s - u ■ y + <p(v)dy 

for (s, w, z) G R x S"- 1 x [0, (2M)" 1 ]. If we put p = s - u ■ y + \y\ 2 z/2, then 
V y p = -co + zy. Since \V y p\ > l-z\y\ > 1/2 for z G [0, (2M)" 1 ] and |y| < M, 
Q[<p] can be written as an integral of a compactly supported function ip over 
a hyper-surface {y G R n ; s — u ■ y + \y\ 2 z/2 = 0} which smoothly depends on 
(s, w, z) G R x S"- 1 x [0, (2M)" 1 ], and we see that 

g[(p) g c°° (r x § n - x x [0, (2M)- 1 ]) . 

We also see that if G[tp](s, u,z)^0 for some (w, 2) G S™" 1 x [0, (2M)" 1 ], then 
we have — 5M/4 < s < M. Indeed, since the assumption implies s = u ■ y — 
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(\y\ 2 /2)z for some y and (u,z) with \y\ < M and (u,z) G S™" 1 x [0, (2M)" 1 ], 
we get 

5M . \y\ 2 1 lyl 2 . . . 

— ~\y\ — ^7777 — u ' y — 7r z \ = s ) — \y\ — m. 



4 ~ 2 2M 

For a multi- index a with lotl = k < 1, we have 



G[dZ<p](s, u, z) = d k s 5^s - cu ■ y + K 2 j (w - zyT^y)dy. (3.10) 



We define 



% w [^]((7,u;,z) := 2 ( 27r )(m-i)/2 J^*+ ^ ( s - (7 + y^J 0M(s,w,*)ds. 
Then we obtain 

r<- m -W(E m {t)*<p){x) = n m [<p](<r,u,r- 1 ). (3.11) 

Since we have 



n m [(f](a,u,z) 
("l) fc 



X+ +(1 " m)/2 ^ - ° + ^-z) (8tQ[<p\) (s, w, z)ds 



2(27r)( m - 1 )/ 2 7 R A+ V 2 

for any nonnegative integer k, and since we have x+ G C 1 (M) for a > 1, we 
can easily see that Um[<p] eC°°(Rx S"" 1 x [0, (2M)" 1 ]). Moreover we have 

a < M in supp"H m [<^] when m is even, (3-12) 

and 

\a\ < M in supp"H m [(^] when m(> 3) is odd. (3.13) 

Indeed <KTM and (T3TT3]) for 2 ^ follow from (ET5)) . Q) , and (13TT]) . while 
they follow immediately from the definition of H m when 2 = 0. 

3.2. The Radon transform and the Friedlander radiation field. Let 

iS(lR n ) denote the set of rapidly decreasing functions on M. n . For ip G iS(lR n ) 
we define the Radon transform lZ[<p] of ip by 

K[<p](a, cj):= f <p(y)dS(y), (^)GRxr 1 , (3.14) 

where dS(y) denotes the surface element on the hyperplane {y G M n ; y ■ u = 
cr}. It is easy to see that TZ[<p] G <S(R x S n_1 ). Since we have 1Z[(p](a,u) = 
G[(p](a,U),0), it follows from f )3.10p that 

K[d£<p](a,u>) = u a d k a K{v](a,u) (3.15) 
for any multi- index a with \a\ = k < 1. 
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For a positive integer m and <p G 5(IR n ), we put 
TZ m [<p] (ex, u) : = 2 ( 27r ^ m -i)/ 2 f R X+~ m)/2 (s - o)K[(p] (s, uo)ds 

= 2(27r) Li)/2 (X-" mV2 * ^MO, G R x S- 1 , 

where := cr ) f° r a G C, and * is the convolution with respect to 

(x-variable. Note that we have lZ m [<p] = U m (-, -,0) G x S™" 1 ). For any 

multi-index a with |a| = k < 1, we obtain from ( 13.15P that 

1Z m [d^](a,u) = 2(27r) e m -i)/2 jj±- m)/ \s - a)u«d k s n<P\Mds 

( — l) k U a f (l- m )/2-fc/ s^r w \ 7 

= 2( 27 r)(-^ / K X+ (* 

=u^ft m MM- (3.16) 
If ^(x) — for |x| > M, then we immediately see by (13. 14ft that 

K[(p](<r,u)) = 0, \a\ > M, u G S" _1 . (3.17) 
Consequently we get 

K m [cp){a,co) = 0, cr > M, u G § n_1 . (3.18) 
When m(> 3) is odd, since suppx+ m ^ 2 = {0}, we obtain 

ll m [ip}{(j,uj) = 0, \a\ >M, co G § n-1 . (3.19) 

Lemma 3.1. Let m be a positive integer with m > 2. For (p G C^°(lR n ) ; a 
nonnegative integer j , and a multi-index a with \a\ = k < 1, there is a positive 
constant C such that 

\diTZ m [d^]{a,u)\ < C (a)- j - k+(1 ~ m)/2 , (a,u) Glx S"" 1 . 

Proof. Suppose that ip(x) =0 for \x\ > M with a positive constant M. Since 
diK m [d%<p) G C°°(R x S"- 1 ), we get 

|«9^ m [«9>](cr,u;)| < C < C (2M) J+fc+(m - 1)/2 ( (7 )-^-fe+a-H/2 (3>20) 

for (cr, w) G [— 2M, M] x § n_1 , which implies the desired result for odd m 
because of ( I3.19p . 

Let m be even. In view of ( 13 . 18[) and (I3.20p . it suffices to consider the case 
where a < —2M. Then we have s — o > \o~\/2 > C (cr) > for \s\ < M. Hence 
we obtain from ( 13. 16[) that 



\diK m [d«Lp]{a,u) 



( i \ j+k, ,a r M 

{-l) J CO / -j-k+(l-m)/2 ( \<r>\ M u 

2(27r)("-D/* j_ M X+ (s-a)n[ V }(s,u)ds 

/M 
(s - a)-^ k+ ^- m ^ 2 ds < C ( a )-i- k +(i-™)/2 ^ 
-M 



because of (I3.17p . This completes the proof. □ 
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Lemma 3.2. Let n > 2, and m be a positive integer with m > 2. Suppose that 
ip G C^°(]R n ) and (p(x) = for \x\ > M with a positive constant M. Then, for 
any integer j with < j < 2, and a multi-index a with \a\ = k < 1, there is a 
positive constant C such that 



^l 2 dl{E m {t) * dy)(x) - ((-^7U5>])(r - t, 



CO 



< Or' 1 (t - r)- ] - k+{ *- m )/2 , r>->2M, (3.21) 
where r = \x\ and u = \x\ x. 

Proof. Recall the definitions of Q and H m in the previous subsection, and that 
we have 

g[d^}( s ,u,o) =n[d^]( s ,u), n m [d^}(a,u,o) = n m [d^}(a,oj). 

We suppose that r > t/2 > 2M. By (I3.8P and (I3.18p . we may also assume 
r <t + M. 

First we assume that j = 0. We put a = r — t as before. By (13. lip we obtain 

\r^' 2 (E m (t) * - (K m [d a M){r - t,u)\ 

= iHmid^^^r- 1 ) - n m [d:<p}(a,u,0)\ 

<r~ l ! \{d z U m [d a x v\)(°^,9r- l )\d9, 
Jo 

which leads to (I3.2ip with j = if we can show 

\d z U m [d^}(a,u,z)\ < C (a)- kH3 - m)/2 (3.22) 

for {a,u,z) G lx S"- 1 x [0, (2M)" 1 ]. Since U m {d%y] G C°°(R x S"" 1 x 
[0, {2My 1 }), we have 

\d z n m [d:^(a,u,z)\ < C < C (a)- k+(s - m)/2 (3.23) 

for (<r,u,z) G [-2M,M] x S"" 1 x [0, (2M) -1 ], which leads to (13T22|) for odd 
m because of (I3.13p . Let m be even. In view of (I3.12p and (I3.23p . it suffices 
to show (I3T22|) for (a,u,z) G (-oo, -2M] x S"" 1 x [0, (2M) -1 ]. Suppose a < 
-2M, u G S"" 1 , and < z < (2M)" 1 . We compute 

2(2nY m ~^ 2 d z n m [d:if](a,u,z) 

= I y x; (1+m)/2 (s-a+ £s) ( a , w , 

+ jf X?^ )/2 " or + y ^ (d z g[^))(s, CO, z)ds =: /j + J 2 . 



I'll < 



$ a (s,u,z) 
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Since s - a + a 2 z/2 > -5M/4 + \a\ > 3\a\/8 > C (a) > for a < -2M, 
z > 0, and s > — 5M/4, recalling (I3.10p we obtain 

where 

J 5 (^s-u-y + (w - zy) a Lp{y)dy 

which is a C°°-function in R x S^ 1 x [0, (2M)" 1 ] with -5M/4 < s < M in 
supp$ Q . Similarly we get I-Z2I < C (o-) _fc+ ^ 3_m ^ 2 . This completes the proof of 
(13.221) for even m, and (13.211) with j = is established. 
Note that by (13.211) with j = and Lemma 13.11 we get 

y m -^l 2 (E m (t) * d» (x) I < C ( (T )- fc +( 1 - m )/ 2 . (3.24) 

Direct calculations lead to 

d t E m (t,x) = 2<KtE m+2 (t,x), -d a K m [d*ip\{a,uj) = 27rTZ m+2 [d" ip](a,uj). 

Hence it follows from (13T24D and fl3~2TD with j = that 
| r(m -i)/ 2(9i ^ m(t) + d> ^ %) _ (-d^K^dy]^^ 

< 2n\r^/ 2 {E m+2 (t) * «9»(x) - ^ m+2 [9>](a,a;)| 

+ 2ny m -^ 2 (t - r)(E m+2 (t) * <9»(x)| < Cr" 1 (j)"^ , 
which is (13.211) for j = 1. Observing that 

<9 2 £ m (t,x) = 27rE m+2 (t 7 x) + (27rt) 2 E m+A (t) 

and 

(-9 CT ) 2 ^ m [^>](a,a;) = (27r) 2 ft m+4 [<9>](a, w), 
we can show (13.211) for j = 2 in a similar way. □ 

For (y?, t/>) G (5(R n )) 2 , we define the Friedlander radiation field J~o[(f, ip] by 
Jo[^,^](^w) = -^ n [^]((T,a;) +72 n [V'](«7,a;), (a,cu) eRx S 71 " 1 . (3.25) 
Observe that this definition is a generalization of the previous definition (12.31) 
that was given only for n = 3 and (<p,ip) G (C^°(IR 3 )) , because TZ 3 [h] = 
(47r) _1 7?.[/i]. The next lemma is a slight refinement of [61 Theorem 6.2.1] (see 
also Friedlander [3] and Katayama-Kubo |12j). 

Lemma 3.3. Letn > 2 anrfc > 0. Letw be the solution to the Cauchy problem 
(13.21) - (13.31) . If (w , u>i) G (C£°(R n )) , i/jen i/jere zs a positive constant C such 
that 

\r^ 1)/2 w(t, x) - W(r -ct,u)\<C(t + r)~ l (ct - r f~ n)/2 (3.26) 

and 

\ r (n-i)/2 dw ^ x j _ Co c ( x )(d a W)(r - c t,u)\<C{t + r)- 1 (ct - r ) {1 - n)/2 (3.27) 
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for all (t, x) G [0, oo) x (R n \ {0}), where W(a, u) = ^[wo, c~ lr wi](a, u), oo c (x) 
is given by (12.21) . r = \x\, u = |x| _1 a; ; and d = (do, di, . . . , d n ). 

Proof. We may assume that c = 1, because the general result is easily obtained 
by a change of variables. We assume that Wq(x) = Wi(x) = for \x\ > M. 

Firstly we suppose that r > t/2 > 2M. Then we have r _1 < C (t + r}~ . Let 
| a | = k, and suppose that j and k are nonnegative integers with < j + k < 1. 
From (13.41) we get 

#a>(t, x) = &> +l (E n (t) * <9> ) (x) + d{ (E n (t) * 9>i) (x). 
Then Lemma 13.21 (with m = n) and (13. 16[) lead to 

r {n - 1)/2 dld«w(t,x) =(-l) j co a (di +k W)(r -t,u) 

+ 0({t + r)- 1 (t - r )-i- fc +(3-«)/2 j ; 

which implies f!3.26|) when j + k = 0, and (I3.27P when j + k = 1. 

Secondly we suppose that either r < t/2 ot t < AM holds. Then we have 
(t — r) _1 < C (t + r)^ 1 . Hence by Lemma [3.11 (with m = n) we get 

|(fl£W0(r-*,a;)| < C (t - r}- jH1 - n)/2 < C{t + r}- j+(1 - n)/2 . (3.28) 

Now we are going to prove 

\A n -^ 2 d a w(t,x)\ <C(t + r}- lal+{1 - n)/2 (3.29) 

for (t, x) with either r < t/2 or t < 4M, which implies (13.261) and (I3.27P with 
the help of (I3.28p . If we assume — 2M < r — t < M in addition, then we get 
t < 4M and r < 5M. Now, from (13.41) we can easily obtain 

\ r <n~x)/2gp w fo x )\<c <C(t + r )-l Q l+( 1 -")/ 2 ; t < 4M, r < 5M. 

When n is odd, this shows (13.291) because of (13. 5p and (13.61) . Thus we assume n 
is even and r — t< —2M. Accordingly we have t/2 > r, and (13.71) immediately 
implies the desired result. This completes the proof. □ 

3.3. The translation representation. For cp G iS(R n ), we write (p(— J r [ l ^\) 
for its Fourier transform. To be more precise, we put 

m = •%>](£) = 7^r 2 [ e**ip{x)dx, £ G R B , 

where i = y/—l. For a function ip = ifj(o~,u) G <S(R x S n_1 ), we define 
~ ~ If 00 

i){p,u) = J r [i)]{p,u) := -= e- lpa ip{a,u)da, (p, u) G R x S n_1 , 

V27T J-oo 
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which is the one-dimensional Fourier transform of ip(-,u) with a parameter 
uj G § n_1 . Let <p G S(R n ). For (p,u) eRx $ n -\ we have 

1 —n&}( p ,u). (3.30) 



(2tt)(«- 1 )/ 
Equation ( 13 .30 j) implies 

K&\(-p, -u) = (2vr) (n - 1) ,/2 fi(puj) = K\y\(p, oo) } (p, u) G R x S^ 1 . (3.31) 

In other words lZ[tp] is an even function in (p,u). 

Lemma 3.4. For ((p,ip) G (iS(IR n )) 2 we /icrae 

||d ff J r o[y,V , ]|U a CKxS«-i) = ll(^^)lko(ffi")- 

Proo/. Since we have K n [h\(a, uj) = B n ( X { l~ n)/2 *ll[h}(-,uj))(cr) for h G <S(M n ) 
with B n = 1/ (2(27r)( n - 1 )/ 2 ), we get 

K^h](p,u>) = F x [n n [h){;u)]{p) = V^B^x^^^nh^u), (3.32) 

where T\ denotes the one-dimensional Fourier transformation (of a tempered 
distribution). It is known that 

Jkb£](p) = -Le^+^p+^or- 1 

V Z7T 



for every a G C, where z fe for z G C \ (— oo,0] and 6 G C is given by z b = 
exp(61og2;) if we write logz = log \z\ +iaxgz with — 7r < argz < tt, and 
(p + i0) b is defined by (p + i0) b = \im £ ^ +0 (p + ie) b (see Hormander 0, Example 
7.1.17] for instance). Especially we have 

\Mx ( -' n)/2 ](p)\ 2 -^\pr 3 , peM. 

Therefore, for <p, ip G iS(M n ) and nonnegative integers j, fc, it follows from the 
Plancherel formula for the one-dimensional Fourier transform and (I3.32j) that 

(^M,W]) #xrl) 

= {-i) k VBl [ ( T ^]{p,uj)^]{p,uj)\p\ n -^ +k dp\ dS^, (3.33) 

where dS w denotes the surface element on § n_1 . By (I3.3ip . we see that the in- 
tegrand on the right-hand side of f!3.33|) is an odd (resp. even) function in (p, uj) 
if j + k is odd (resp. even). Hence we have (d%lZ n [ip], do-lZnty]) L2 r Mx §n-i\ = 0. 
Now we find 

\\daJ 7 a[<P, V , ]l|i2( RxS n-l) = ||^^7lM||i2( KxS n-X) + 1 1 d a 1Z n [ip] 1 1 £2( Rx gr»-l) • 
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We obtain from f[3733|) . f l33Tj) . and f l33U|) that 

2 p"" 1 ^ j dS u 



ll«nMlli 2 (Rx^- 1 )=2^ / / 



1..-M.2 I 



Just in the same manner, we obtain 



This completes the proof. □ 



By Lemma 13. 4[ we can uniquely extend the linear mapping 
(<S(R n )) 2 3 fai/)) ^ d a T [cp,i)} G L 2 (R x S"- 1 ) 
to the linear mapping T from Ho(M. n ) to L 2 (R") with 

||T[^V]||^(Rxs»-i) = \\M)\\h <b»), M) E H (R n )- (3-34) 

This mapping T is called the translation representation in Lax-Phillips [20J. 

The following lemma was essentially proved in [20] for odd n, and in [19] for 
even n (see also Melrose 



Lemma 3.5. The mapping T defined above is an isometric isomorphism from 
H (R n ) to L 2 (R x S"- 1 ). ' 

Proof. What is left to show is that T is surjective. Because of (I3.34j) . we 
only have to prove the following: For any function v in some dense subset of 
L 2 (R x S n ~ l ), there is (<p,ip) G H (W n ) such that T[<p,ip] = v. 

We write h G 5 (Rx§ rt_1 ) if h G 5(Rx§™ _:L ) and there is a positive constant 
5 such that h(p,u) = for all (p,u) G {-5,5) x S™" 1 . We put 

5i(M x S"- 1 ) = \?- x [h\] h g S (R x S"- 1 )} , 

where 

~ If 00 
F-\h](o-,u) = -= / e^h(p,u)dp, 

which is the one- dimensional inverse Fourier transform of h(-,u) with a pa- 
rameter io G S n_1 . It is easy to see that S (R x S n_1 ) is dense in L 2 (R x § n_1 ), 
and hence S X (R x S™" 1 ) is also dense in L 2 (R x S"- 1 ). 

Let v G «Si(R x S"- 1 ). We want to find (<p, ip) G H (R n ) such that T[<^, V] = 
v. We put 



«o(P, W) = - r (1 _ n)/2l ; W), (p, W) e R X 



5n— 1 



2^B n ^[ X ( l- n),1 ]{p) 
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Note that we have v G S (R x S™" 1 ), because v G S (R x S™- 1 ), and the 
singularity of l/J-i [yL 1 ™^ 2 ] lies only at p — 0. For (Gi n \ {0} we put 



We also set vi(0) = v 2 (0) = 0. Then, since v G 5 (K x S"" 1 ), we find that 
t>i,t>2 G 5(M n ). Hence if we set <p = •7 r ~ 1 [fi] and if) = T 7-1 ^], then we get 
(</?,-0) G (5(M n )) 2 C H (R n ). Using f l33U|) and f l3T32j) . we get 

7 [d^o [(f,if)]\ (p, u) = p 2 TZ n [if] (p,uj) + ipTZ n [if)] (p, u) = v{p, u) 

for (p,u) G K. x S n_1 , which shows T[(f,if>] = daF^^] = v. This completes 
the proof. □ 

Theorem I2.ll is an immediate consequence of the following lemma. 

Lemma 3.6. Let c > 0, and (wo,Wx) G H (R n ). Let w be a solution to the 
Cauchy problem (I3.2l) -( !3~3l) . and 



W(t,x) = \x\- {n - 1)/2 T[w , c^wijdxl - ct, Ix^x) 
for (t,x) G [0,oo) x (R n \ {0}). Then we have 

^||^Mt,0-^(OW(t, 011^ = 0, 

where := (c _1 <9 t , <9i, . . . , d n ) and cu(x) := toi{x) = (—1, |a;| _1 a;). 
Proof. Let e > 0. There is {w* ,wl) G (C °°(M n )) 2 such that 



Let w* and W* be defined similarly to w and W, respectively, by replacing 
(wq,Wx) with (wq,w*). Recalling (|3.l|) . we obtain from the energy identity 
that 



«2(0 



«i(0 



2( 27r )(n-l)/2^|2 

2i(27r)(™- 1 V 2 |£| 




i||^(t,-)-^V(t,-)||i 
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Since \fi(x)\ = 2, using (I3.34j) we get 

lp(-)(mr)-mtr))\\l {wn) 

\T[w - Wq,c~ 1 (w 1 - wl)](r - ct,u)\ 2 dS u ^j dr 
< lirK-w^c-^Wi-w^Hlla^sn-i) 

= ||(w ,C~V) - K,C~X)lltf (M") < 

By Lemma [3.31 we obtain 
L ~' c) w*(t,.)-u(-)W*(t r )\\l 2m 



2 



2 

< 

2 



< \ I I I \r in - 1)/2 d ic) w*(t,ru) - uj(ru)(d a ^)(r - ct,oj)\ 2 dS„ ) dr 



V./S 

oo 



poo 

<C (1 + t + r)" 2 (l + \t- r\)-^- 1] dr < (7(1 + t) -1 , 

where J-"q (a, w) = J-o[wq, c~ 1 u;5 s ](cr, a;). Hence there is a positive constant t > 
such that t >t implies 

^ii^V(t,-)-^(-)w*(t,-)iiV)< (|) 2 . 

To sum up, for any e > there is a positive constant £ such that we have 

-L\\d^w(t,-)-u(-)W(t,-)\\ L2{Rn) <£, t>t . 

This completes the proof. □ 

Now we are in a position to prove Theorem 12.11 
Proof of Theorem 12.11 Suppose that there is V G L 2 (M x § n_1 ) such that 

lim||^(t, •) - uj c (-)V% -)\\ L2m = 0, (3.35) 

where 

V\t,x) = \x\- {n - l),2 V(\x\ - ct, (3.36) 

Then we get 

\\m\\d^v(t,-)-^(-)VKt,-)\\ LHRn) = 0, (3.37) 

where and cU(:r) are defined as in Lemma 13.61 We define (wq, c~ 1 wi) = 
T _1 [^](e H (R n )), and let w be the solution to the Cauchy problem (1X2]) - 
(13.31) for this (wo,^i). Then it follows from Lemma [3.61 that 

lim \\d^w(t, ") - <3{')V*(t, -)\\ L2{Rn) = 0, (3.38) 



which, together with ( 13.371) implies 

lim \\v(t,')-w(t,')\\ E ,c= \im ^=\\d^v(t r ) - d^w(t,-)\\ L 2 {Rn) = 0. (3.39) 
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Conversely, suppose that there is (w ,Wi) G H {R n ) such that f!3.39|) holds 
for the solution w to the Cauchy problem (I3.2p -f l3"!?j) . If we put 

V(a,u) = T[w Q ,c- l Wl }{a,uj){e L 2 (R x S"" 1 )) 

and define V$ by ll£5B|) . then Lemma EH implies fEQ8|) . Now ( l3~39l) and 
yield (13.371) . which is equivalent to (13. 35ft . This completes the proof. □ 

4. Vector Fields Associated with the Wave Equations 

We restrict our attention to the three space dimensional case from now on. 
We introduce vector fields 

S -=td t + x-V x , 

Q =(Qi, Q 2 , ^3) := x x V x = (x 2 d 3 - x 3 d 2 , x 3 d x - xxd 3 , x x d 2 - x 2 dx), 

where the symbols • and x denote the inner and exterior products in M. 3 , 
respectively. We also use d = (d t , V x ) = (<9 a ) < a <3- We define 

z = (z , Zi,..., z 7 ) = (s, n, d) = (s, n l5 n 2 , n 3} d 0} d u d 2 , d 3 ). 

These vector fields are compatible with the system of wave equations having 
multiple propagation speeds, since we have [□<., flj] = [D c , d a ] = for 1 < j < 3 
and < a < 3, and [D c , S] = 2D C , where c > 0, and [P, Q] = PQ - QP 
for operators P and Q. Using a multi-index a = (a ,ai, . . . ,07), we write 
Z a = Zq° Z^ 1 ■ ■ ■ Zj 7 . For a nonnegative integer s and a (scalar- or vector- 
valued) smooth function <p = ip(t,x), we define 

\<p(t,x)\.= \zMt,x)\. 

\a\<s 

We can check that we have [Z a ,Z b ] = ELoQ' Z d and [Z a ,d b ] = YL Q Dfd d 
with appropriate constants and D^ b . Hence, for any multi-indices a and 
(3, and any nonnegative integer s, there exist some positive constants C a ^ and 
C s such that we have 

\z a z^(t,x)\ <C a ,p\<p(t,x)\ la]+m , 

3 

C; l \d V {t,x)\ s < ^^\d a Z a <p(t,z)\ <C s \d V (t,x)\ s 

\a\<s a=0 

for any smooth function ip. 

We write r = \x\, u = (ui,u 2 ,u 3 ) = and d r = Y^j=i u 3®y Then we 

can write S = td t + rd r . For c > 0, we define 

«9i c) =d t ± cd r , and l4 c) = ±^-d£ c) = l -(d r ± 0-%). 

Since we have d t = -c(D ( f> - Df), d r = + D+\ and 
(1 + r)d ( + ] =cS- (ct - r)d t + (d t + cd r ) } 
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there exists a positive constant C such that 

\(d t - {-c)D®)<p{t,x)\ + \{d r - D®)<p(t,x)\ 

< {\Z<p(t,x)\ + \ct - r\ \d<p(t,x)\) (4.1) 

1 + r 

for any smooth function ip. Since we have V x = u>d r — r _1 w x Q, we get 

\(d k -u k d r )<p(t,x)\ < C{l + r)- l \Ztp{t,x)l k = 1,2, 3. (4.2) 
From (14. ip and (14.21) we obtain 

C 

dcp(t, x) — u c (x )D ( L ) (p(t,x) \ < (\Ztp(t,x) \ + \ct-r\ \8<p(t,x)\) (4.3) 

1 + r 

for any smooth function ip, where u c {x) is given by ( 12.2ft . As an immediate 
consequence we also get 

rd(f(t, x) — oo c (x )D (c) (r<p(t, x))\<C (\<p(t, + \ct - r\ \d<p(t, x)\) . (4.4) 

We remark that the term \ct — r\ \d(p\ was not needed in the estimates used in 
[ID] instead of ( 14. 3 p and ( 14.4ft : however the vector field L = tW x + xdt, which is 
not compatible with the multiple speed case, was involved. We need the term 
\ct — r\ \d(p\ here to compensate the lack of the vector field L. These kinds of 
identities and estimates without the vector field L were developed and used in 
PU, CP, [23], [21], [29], and so on (see also [9J, p], [21], [25], and [26] for the 
related topics). 

Using (14. 3p . we can easily show the following estimate for the null forms Qq 
and Qab given by (11. 7ft and (11. 8ft (see [23], [27] and [29] for the details of the 
proof): 

Lemma 4.1. Let c > 0. T/ien we /iave 

|g (^,^c)|+ \QobM)\ <c (r)- 1 (\zcp\ \fri>\ + |^||^|) 



0<a<b<3 



+ C (r)' 1 (ct - r) \d<p\ 



at (t, x) G (0, oo) x E 3 with r = \x\. 

Outline of proof . Fix c > 0, and we define R = (R a )o<a<3 = d — u c (x)D ( f\ 
Substituting d = uj c {x)D_ + R, we get 

\QoM; c)| + ^ |Q„ 6 (^, V)| <C(|£>1 C V| W + 1^1 |^V| + \R<p\ W\) 

a,b 

<C (\d<p\ + |%| |^|), 

where we have used l-R^I < to obtain the last line. The point here is 

that terms including (£)^<£>) (D_ "0) are canceled out because of the structure 
of the null forms. Now, using (14.31) to evaluate \R(p\ and \Ri[>\, we obtain the 
desired result. □ 
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Let ci, . . . , cjy satisfy (jl.5p . We define 

A =|(t,r)eR 4 -xR + ;r>^>l|, (4.5) 

where M + = [0, oo). For 1 < j < N, we define 

toj(o-) := max{-2(2 Cj - Ci) _1 cr, 2c7/} (4.6) 

so that we have Uo-grU^ c jt + * — ^oj(°~)} = A . We also define 7~oj(er) := 
Cjioj (°~) + °~ f° r o G M. The following lemma is a modification of a key lemma 
in [10] to obtain the asymptotic behavior. 

Lemma 4.2. Assume (jl.5p . and G {1, . . . , iV}. Suppose the following: 

• fij > 1 and > 0. 

• /Xfc > and Kfc > 1 /or any = 1, . . . , N with k ^ j. 
If v = v(t,r,u) G C l (A x § 2 ) satisfies 

d^ j) v(t, r, w) = G(t, r, w), (*, r) G A , oo G § 2 , (4.7) 

and G satisfies 

N 

\G(t,r,uj)\<Y,B k (t + r)-^ {c k t-r)~ K \ (t,r) G A , w G § 2 (4.8) 
fe=i 

wt/i some positive constants B\, . . . , Bn, then there exists a positive constant 
C such that 

\v(t, r, uj) - V(r - Cj t, u)\ <CB 3 (t + r}' H+1 ( Cj t - r)~ Kj 

+ C B k (t + r)-^ (4.9) 

l<k<N 

for any (t, r) G A and u 6 § 2 , where V is defined by 

POO 

V(a,u) = v(toj(a),r 0t j(a),u) + / G(s, CjS + a, oo)ds. (4-10) 

J t j (a) 

The constant C above is determined only by c k , Pk, and n k with 1 < k < N. 
Proof. First we note that V is well-defined because of ( 14. 8p . By (14 .7p we find 

v(t, Cjt + a, oo) = v(t j(a), r j(cr), oo) + / G(s, CjS + a, uS)ds, t > t j(o~) 

Jto,j(cr) 

for (a, oo) G R x § 2 . From ( OU|) and (gSJ we get 

/.oo ^ 

-V(a,u)\< / |G(s, Ci s+o-,a;)|ds < Cj2 B kh(t,a), (4.11) 

where 

^ *)= I (1 + r + P)^(l + |c fc r - p\)-« \ {TtpMsiCjS+a) ds 



ASYMPTOTIC BEHAVIOR FOR SYSTEMS OF NONLINEAR WAVE EQUATIONS 21 

for 1 < k < N. By direct calculations, we obtain 

«*■ " ] = (c, + l)(„,-l) (1 + {C ' + ^ + a) ~" +1(1 + kl) " % <4 - 12) 
since fij > 1. If k ^ j, then we have 

/CO 
(l + \(c k - Cj )s-a\)- Kk ds 



OO 



<(1 + (C, + l)t + or)-^ 1 / (1 + \T\)- K HT 

Cj\ J—oo 

£ (^lk^ (1 + fe + 1)t + ^" (4 ' 13) 

for t > t j((r) and cr G M, because we have > 0, > 1, and c& 7^ Cj. 
Now we obtain (14. 9p immediately by putting a = r — Cjt in ( 14. lip . (I4.12p . and 
(14131) . □ 

5. Proof of Theorem 12.21 and Corollary 12.31 

In the following, we write r = \x\ and u = Let the assumptions in 

Theorem 12.21 be fulfilled. Then we have the global solution u to the Cauchy 
problem ( ll.ip - fll.2p for sufficiently small e by the global existence theorems in 
[27], [28], or [29]. For j = 1, . . . , N, we put 

Uj(t, x) = Uj(t, x) — evPAt, x), 

where w° is the solution to □ Cj w° = with initial data w° = fj and dtu® = g-j 
at t = 0. From flTT]) and fO]) . we have 

□ CjM j = ^(dw, <9 2 w) in (0, oo) x M 3 , (5.1) 

u)$,x) = (dtUj)(0,x) = 0, x G M 3 . (5.2) 

We put 

Suppose that we have f(x) = g(x) = for \x\ > M. Then, as in (I3.5p . ( 13. 6p . 
and (13.191) . we get 

u °.(t, x) =0, \r - Cjt\ > M, (5.3) 

Uj(<r,u) =0, |a|>M (5.4) 
for 1 < j < N. Hence Lemma 13.31 implies 

\rv^(t,x) - £/°(r - c,-t,x)| + |rdu°(t,x) - u; Ci (x)(9 CT C/?)(r - Cjt,x)\ 

< C (t + ry 1 ( Cj t - r)' 1 , l<j<iV (5.5) 

for G [0, oo) x (R 3 \ {0}), because (JED and (JED imply 1 < (c^ - r) < 
(M) on the support of the functions on the left-hand side of (15.51) . 
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Let A be the set of (t, x) G [0, oo) x M 3 with (t, r) G A , where A is given 
by (j4.5P ; namely we put 



(5.6) 



A= G [0,oo) x (M 3 \{0});|x| > ^ > 1 

We set A c = ([0, oo) x (R 3 \ {0})) \ A. Recall the definitions of d ( ± ] and D { ± ] 
in the previous section. 

Proof of Theorem \2.2\ (1). Following the proof of the global existence theorem 
in [18] (see also [7J Proposition 4.2] and its proof) we obtain 

\u k (t, x)\ 2 + e- l \ul(t, x) | 2 <Ce (t + r)" 1 log (l + l±M±Il ) 

|0M*,*)li + e-^^^x)!! <Ce (r)- 1 ( Cfc * - r)" 1 (5.7) 

for 1 < A; < A/". We choose small 5 > 0. Then, by (15. 6 p we get 

\u k (t,x)\ 2 + e~ l \u\{t,x)\ 2 <Ce(t + r)- 1+s (c k t - r)~ 5 (5.8) 

for 1 < k < N, because there is a positive constant C§ such that we have 
log(l + z) < C s z 5 for z > 1. 

Fix j = 1, 2, . . . , N. Switching to the polar coordinates, we get 

r{D c .u}){t, rco) = d { + j) d { ° j) [ru)(t, ru)) - cJr^A^t, ru) (5.9) 

for (t, r, u) G (0, oo) x (0, oo) x§ 2 , where A w = Y^k=x * s the Laplace-Beltrami 
operator on § 2 . We put 

Vj{t,r,uj) := Di^frujforw)). (5.10) 
Then and (JEl lead to 

d { + j) Vj{t,r,u) =Gj(t,r,u}), (5.11) 

where 

Gj(t,r,u)) := -— {c^A^ftrw) + rF j (du(t,ru),d 2 u(t,ru))} . 

We suppose (t, r) G Ao and u£§ 2 for a while. Note that we have 

(1 + 2c^ i y 1 (t + r) <r < (r) < (t + r) , (t, r) G A . 

From (I5.8p we obtain 

\r- l ^u){t,ruj)\ < Ce 2 (t + r}- 2+5 { Cj t - r}- 5 . (5.12) 

Recall the definitions of N j} R], and Rf given in (ODD . ffLTTj) . and (Q2jl . 
respectively. By Lemma |4~T| (15. 7p . and (15. 8p . we get 

r|A^(«9w,<9 2 «)| <Ce 2 (t + r}- 2+5 {cjt-r}- 1 - 5 . (5.13) 

Noting that we have 

(c k t - r) (eft -r) >C (t + r) min{(c fc t - r) , (c t t - r)} 
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for Cfc 7^ q, we obtain from (15. 7p that 

r\R]{du,d 2 u)\ <Ce 2 (t + r)~ l ^ (c k t - r)~ l (c t t - r)' 1 

k^l 
N 

<Ce 2 (t + r)- 2 ^{c k t-r)- 1 . (5.14) 
fc=i 

By (15.71) we also have 

r\Rf(du,d 2 u)\ < Ce 2 J2{t + r )~ 1 (c k t - r)~ 2 . (5.15) 

Now fl5TT2l) . fl5TT3|) . (Oijl . and fl5TT5|l lead to 

|G,(t, r, w)| <CV 2 (t + r>" 2+5 ( Cj t - r>~ 5 

+ CV 2 ^ (t + r)" 1 ( Cjfc t - r>^ 2 . (5.16) 

We define 

/■oo 

V^-(cr, w) = Uj(to,i(o')>7'Q,i(c)j w ) + / Gj(s, CjS + cr, u)ds, (5.17) 

J t ,j (ct) 

where tojC ") is defined by (I4.6p . and r j(a) = Cjt j(a) + a. Since we have 
(15. lip and (I5.16p . Lemma [4.21 implies 

\ Vj (t, r, u) - Vj(r - Cj t, u)\ <Ce 2 ((t + r y 1+S ( Cj t -r)' 5 + (t + r)" 1 ) 

<Ce 2 (t + r)- 1+s ( Cj t-r}- s (5.18) 
for (t,r) G A and u G § 2 . By f l5U0|) . (l5T7jl . and (Q|l . we get 
\vj(t,r,u)\ <C (r\du){t,ru)\ + \u){t,ruj)\) 

<Ce 2 (( Cj t - ry 1 + (t + r}- 1+S ( Cj t - r>~ 5 ) 

<Ce 2 ( Cj t - ry 1 . (5.19) 
Hence, putting r = Cjt + a in (I5.18p . we get 

\Vj(ff,u)\ <\ Vj (t, Cjt + <r,u)\ + CV 2 (( Cj + l)t + a}^ 1+<5 (a)" 5 

^((a)- 1 + (( Cj + l)t + a)- 1+5 (a)- S ) (5.20) 

for all (a, u) G R x S^ 1 and i > toj(o")- Taking the limit in (ODD as t ^ oo, 
we obtain 

\V j (a,u)\<Ce 2 ((r)~ 1 , (ff,w)Glx§ 2 . (5.21) 

Now we define Pj(a,cu) = e~ 1 Vj(a, u). Then, recalling the definition of Vj, we 
obtain from (JO]), §5J§, and d£2J that 

|r9uj(*,x) - w Cj (x)^(t,r,u;)| < Ce 2 (i + r)~ 1+5 (qt - r) _<5 , (t,x) G A. 
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Therefore (I5.18P leads to 

\rdu){t,x) - £U Cj (x)P}(r - Cjt,u)\ < Ce 2 (t + r)~ 1+<5 ( Cj t - r}~ 5 (5.22) 

for (t, x) G A. Finally we define 

Pj((T,u) = d a U°(a,u) + P}(<r,u). (5.23) 

Then $3} and fl5T22|) yield ([21} for (t, x) G A. Noting that fICT]) and the 
definition of Pj lead to 

\Pj((t,uj)\ < Ce(ay\ (it,w)6Rx§ 2 , (5.24) 

and recalling f !5.4j) . we immediately obtain (12. 5p by (I5.23p . 

Now what is left to prove is ( 12. 4p for (t,x) G A°. Since we have either t < 
2c^ or r < Cit/2(< Cj-t/2) in A c , we get ( Cj t - r) _i < C (t + r)" 1 . Therefore, 
by (15. 7p we get 

\rd Uj (t,x)\ < Ce ( Cj t - ry 1 < Ce (t + r)~ 1+s ( Cj t - r)' 5 (5.25) 
for (t,x) G A c . It follows from (Q|l . (IQi]) . and (l5^3|) that 

< C(ay\ (ff,w)eRx§ 2 , 

which immediately yields 

\Pj(r - Cjt,u)\ < C ( Cj t - r)' 1 < C (t + r}~ 1+5 ( Cj t - r)~ 5 (5.26) 
for (t,x) G A c . From (I5.25P and (I5.26P we obtain 

\rduj(t, x) - eu^Pjir - Cj t, w)| < Ce (t + r)~ 1+s ( Cj t - r)' 5 , (t, x) G A c , 
which is ( 12. 4p for (t, x) G A c . This completes the proof. □ 

Proof of Theorem 12.21 (2). Suppose that R 11 has the null structure, and fix p G 
(1/2, 1). Then following the proof of the global existence theorem in [8] (see 
Proposition 5.2 and its proof in [S] specifically), we have 

\u k (t, x)\ 2 + e- l \u\(t, x)\ 2 <Ce (t + r}' 1 {c k t - r}' 1 , (5.27) 

\du k (t, x)\ 2 + e-^duKt, x)\ 2 <Ce (r)" 1 (c k t - r)" 1- " , (5.28) 

which are better than (15. 6p and (15. 7p . Using decay estimates in [15] (see 
Lemmas 3.2 and 6.1 in [15J), we can further improve (I5.28P in the decay rate 
and we get 

\du k {t,x)\t + e-'ldul^x)], < Ce (r)" 1 (c k t - r)~ 2 . (5.29) 
Indeed, similarly to (I5.13P and (I5.14p . we obtain from (I5.27P and ( I5.28P that 
rlNjU <Ce 2 ((t + r}- 1 (r)" 1 ( Cj t - r)" 2 ^ + (r)" 2 ( Cj t - r y l ~ 2p ) 

<Ce 2 (t + ry 2 w4t,r)- l ~ 2p , 
r\R)\i <Ce 2 (r)" 1 ^ (c k t - r)" 1 "" (c t t - r)- 1 '" 

k^l 

<Ce 2 {t + ry 2 W yt 1 r)- 1 - 2p 



ASYMPTOTIC BEHAVIOR FOR SYSTEMS OF NONLINEAR WAVE EQUATIONS 25 

at (t, x) G (0, oo) x M 3 , where w_(t, r) = min <fc<7v (c^t — r) with c = 0. Here 
we have used the estimates (r) (cjt — r) > C (t + r) u>_(t, r) and 

(r) (c fc t - r) (cjt - r) > C {t + r) 2 w-(t, r), k^l. 

Since R u has the null structure, we see that R n enjoys the same estimate as 
Nj. Now we have proved that 

r|F(^,9 2 M)|! <Ce 2 {t + r)- 2 w.{t,r)- 1 - 2p , (5.30) 

and we obtain (I5.29P from the following estimate, which comes from Lem- 
mas 3.2 and 6.1 in [15] (especially see the estimates (3.7), (6.1) and (6.2) in 
[15]): For c,K,fM> 0, it holds that 

(r) {ct-r) l+K \d<P{t,x)\ 

<C7su P (y) 2+K |^,y)|i| T=0 



+ C sup \y\(T+\y\) x+ *w-{TAy\) x+ %UMT>v)W 

(r,),)e[0,t]xl 3 

Let (t,r) G A and u G S 2 for a while. Recall the proof of (1). We use 
flEZZD and fl5^9|) . instead of (JED (or (Ol) ) and (iO) . Then f l5TT3|) . 
and (15.141) are replaced by 

\r- l /\ w u)(t, ruj)\ <Ce 2 (t + r>~ 2 ( Cj t - r)' 1 , 

r\Nj(du, d 2 u)\ <Ce 2 (t + r)' 2 ( Cj t - r)~ 3 , 

TV 

r\R){du, d 2 u)\ <Ce 2 (t + r)~ 3 - r )~ 2 , 

k=l 

respectively. Since R 11 has the null structure, we can use Lemma I4TT1 to get 
r\Rf(du,d 2 u)\ < Ce 2 J2( t + r r 2 (c k t - r)-\ 

k+j 

instead of (I5.15p . These estimates lead to 

\G^t,r,u)\ < Ce 2 (t + ry 2 Uc 3 t - r)' 1 + ^(c fc t-r>- 3 J (5.31) 

for (t, r) G A and u> G § 2 . Now going similar lines to (15. 17j) through ( I5.24p . 
we can construct Pj satisfying 

\rdu](t,x) - eQ C] (x)Pj(r - Cjt,uj)\ < Ce 2 (t + r)" 1 ( Cj t - r)~ l (5.32) 
for (t, x) G A, and 

\PHa,u))\ < Ce(ay 2 , (a, uj) G R x § 2 . (5.33) 



We define 



/>oo 

Uj(a,u):=- P}(\,w)d\ (5.34) 

J a 
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so that we have d a Uj(a, to) = Pj(cr, uj). Note that the right-hand side of (15.341) 
is finite because of (15. 33ft . and that Uj(a, uj) — > as o — > oo. It follows from 
(15T32D and (153TD that 

\d r (ru){t,ruj)) - e(d a Uj)(r - Cjt, uj) \ < Ce 2 (t + r)' 1 ( Cj t - r}' 1 (5.35) 

for (t,x) G A. Since (I5.27P implies ru^(t,roj) — > as r — >■ oo, from (I5.35j) we 

get 



ru^{t,ru) — eUj(r — Cjt, uj)\ < 



d x (Xu 1 j (t, Xuj) - eUj(X - Cj t,oo))dX 



<Ce 2 J (l + t + Xy^l + lcjt-XD^dX. 

(5.36) 

If r > 3cjt/2, then we have 

POO 

(l + t + A)- 1 (l + |c J t-A|)- 1 dA < C / (l + t + X)- 2 dX < Cil + t + r)- 1 . 

J r 

For r satisfying 1 < C\tj2 < r < 3cjt/2, we get 



/ (l + t + A) -1 (l + | Cj -* — X\)~ l dX <C(l + t)- 1 / (1 + {cjt-XD^dX 

Jr 

POD 

+ c (i + t + xy 2 dx 

J3cjt/2 

<C(l + t + r)- 1 \og(2 + t). 

Now it follows from (I5.36P that 

\ru](t,x) -eUj(r -Cjt,u)\ < Ce 2 (t + r)~ x log(2 + t), (t,x)eA. (5.37) 
For (a, u)g!x § 2 , by (1537]) and (l5T2Tj) we get 

|E#(<t,w)| < Cs((c j + l)t + cry 1 \og(2 + t) + Ce (a)- 1 
for sufficiently large t. Now taking the limit as t — > oo, we obtain 

\Uj(a,u)\ < Ce(a)- X . (5.38) 

Finally, putting 

U j (a,uj) = U°(cj,uj) + U}{<j,uj), 

we obtain (ES) and (EZJ) for (t,x) G A from (153]) . (1Q21 . and (15T371) . (EH} is 
an immediate consequence of (I5.33f) and (I5.38H . 

As in the proof of (1), using (l5"27j) . (15351) . (15331) . an d d53BJ, as well as Q52D , 
we can show that (12. 6p and (12.71) hold also for (t, x) G A c . This completes the 
proof. □ 

Now we are in a position to prove Corollary 12.31 
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Proof of Corollary 12.31 Let the assumptions of Corollary 12.31 be fulfilled. Let 
Pj for 1 < j < N be from Theorem 12.21 (1). Recalling (15 .4p . we see from 
fT23|) that \Pj(cr,uj)\ < C (a)' 1 , which implies Pj G L 2 (R x § 2 ). If we put 
Pj(t,x) = |x| _1 Pj(|x| — Cjt, |x| _1 x), it follows from (12. 4p that 

\\duj(t,-) -eu Cj (-)P?(t,-)\\ 2 L2m 

2, 



o ./:: 



< 



rduj(t, no) — eu Cj (roo)Pj(r — Cjt, u) \ dS^dr 

POD 

Ce 2 / (t + r)- 2+2S (cjt - ry 2S dr < Ce 2 (l + t)' 1 -)• 0, t ->• oo. 



2 
OO 



Now Theorem 12. II (with n = 3,c = Cj,v = Uj, and = ePj) implies that each 
Uj is asymptotically free in the energy sense. This completes the proof. □ 

6. Proof of Theorem 12.41 

Suppose that e is sufficiently small, and let u be the global solution to (12 .9p 
with initial data u = ef and dtu — eg at t — 0. Since the last half of (I2.10p 
follows from (15. 6p for any f,g G C^°(1R 3 ), we will prove the existence of (f,g) 
for which the first half holds. Without loss of generality, we may assume 
A 1 > and A 2 > 0. We may also assume C\ = 1 < c = c 2 . 

We suppose that / = 0, and that g = {g\,g-i) is radially symmetric, i.e., 
gj(x) = gj(\x\) with some function g* for j = 1,2. Then u = (mi,u 2 ) is also 
radially symmetric in x-variable, i.e., Uj(t,x) = u*(t, \x\) with some u*(t,r) for 
j = 1, 2. For r G M, we put 

hj(r) := ^*(|r|), ^(t,r) := ru*(t, \r\), j = 1,2. 

From (12. 9p we obtain 

Vj(t,r) = — hj(X)dX + - I G j (T,X)dx)dr (6.1) 

C i Jr-Cjt C j JO \Jr-Cj(t-T) J 

for j = 1,2, where 

Gi(t,r) =^f{d t u* 2 (t, |r|)) 2 = ^-{d t v 2 (t,r))\ 

G 2 (t,r) =^f(d t ut(t, |r|)) 2 = ^(^(t,r)) 2 . 

Assume that g\ is a nonnegative function. Since we have 
d t v 2 {t, r) —e {h 2 {r + ct) + h 2 (r - ct)} 

+ [ {G 2 (T,r + c{t-T))+G 2 (T,r-c(t-r))}dT, 
Jo 

we get 

d t v 2 {t,r) >eh 2 (r-ct)(>0), (6.2) 
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provided that r — ct > 0. We fix M > 0, and assume that 

C := / \h 2 (\)\ 2 d\ > 0. 
Jo 

We suppose that < a = r — £ < M in what follows. We put 
r-t a _ r + t-M_2t + a-M 

r ° ~ c- 1 ~ c- r Tl ~ c + 1 ~ c+i ' 

so that we have cr = r — £ + r and ctl + M = r + £ — T\. Let 2£ > 
(c — l) _1 (c + 1)M hold, so that we have T\ > r . Then it is easy to see that 
£ C V, where 

V ={(r, A); < r < £, r - £ + r < A < r + £ - r}, 
£ ={(r, A); r < r < n, ct < A < cr + M}. 

Therefore, noting that we have Gi(t,r) > for r > 0, we get 



G 1 (r,X)drdX> // G 1 (r,X)drdX. 
v J Js 

From the definition of G\ and (16. 21) . we obtain 

Gl(T , A)dTdA r c r +M i^zfQiV) dT 



2 J T0 \Jcr A 
,1, 2 /• 1 / /• '' , x , 



Co^xe 2 , en + M C^e 2 , (c - l)(2c£ + M) 
log — > — log 



2c cr + M~ 2c & (c+ l)(2c- 1)M' 

(6.3) 

Now we choose a nonnegative function g\ satisfying f^ +1 hi(X)dX > 1. 
Then we get 

/ + £ii(A)dA > 1 (6.4) 

Jr-t 

for (£, r) with 0<r-£<Mand£>M + l 

From ( 16.11) with j = 1, ( 16. 3D . and ( 16. 41) . we obtain 

> C(£ + e 2 log(2 + £)) 

for < r — £ < M and £ ^> 1. This completes the proof. □ 
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